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Abstract 

>>: 

We use the cross correlation sum introduced recently by H. Kantz to study 
symmetry properties of chaotic attractors. In particular, we apply it to a 
■ system of six coupled nonlinear oscillators which was shown by Kroon et al. 

to have attractors with several different symmetries, and compare our results 
with those obtained by "detectives" in the sense of Golubitsky et al.. 
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1 Introduction 



Spontaneous symmetry breaking is one of the most interesting problems in statistical 
physics. While in the past most work has been done on equilibrium systems, it has 
become clear only recently that very interesting similar phenomena occur also in 
non-equilibrium systems. Indeed, the variety of phenomena observed there is even 
larger than in equilibrium models, even if we restrict ourselves to purely deterministic 
(but in general chaotic) systems as we shall do in the present work. Spontaneous 
symmetry breaking, in particular, can be of two fundamentally different kinds. 
Assume that the equation of motion 



(with x(t) in some manifold M) is equivariant (or "covariant" in the physics liter- 
ature) under some symmetry group G acting on A4, i.e. 



for each element g G G. We assume that F has at least one chaotic or quasiperiodic 



In such a case, orbits x(t) can have the following symmetry properties: 

• (1) They are symmetric under each element G G, i.e. gx(t) = x(t), for each 



• (2) Individual orbits are symmetric only with respect to some subgroup E C 
G, but the attractor and the invariant measure are still symmetric under G. 
Notice that £ can be trivial (having just the identity as unique element), 
in which case the instantaneous symmetry is maximally broken, while the 
symmetry of the attractor is still unbroken. Such a case occurs most naturally 
when the system is composed of two identical sub units and G is just the 
permutation group of two elements interchanging them. The transition from 
item 2 to item 1 manifests itself then in synchronization of the two subunits 
|1|, |2|, |3|. Thus (de-)synchronization is just the simplest version of a sort of 
spontaneous symmetry breaking, and it was stressed in 0] that the breaking 
of other symmetries should follow essentially the same lines. 

• (3) Also the attractor is no longer symmetric under the full group G, but only 
under some subgroup S C G (the "symmetry group" of the attractor). In this 
case one has of course coexistence of several attractors related to each other 
by symmetry transformations g G G/H. Transitions between item 1 and item 
3 are possible |J but less common than transitions either between 1 and 2 or 
between 2 and 3. The latter have been studied intensively in f|, 0, |f. 111 the 
most frequent scenario, the different attractors in the symmetry broken phase 
are disjoint, and symmetry is restored when these attractors touch in one or 
several points. 




(1) 




(2) 



attractor. 



geG. 
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In the present paper we shall only deal with item 3, and with transitions between 
2 and 3. We shall not deal with item 1 or with transitions between 1 and 2. More 



precisely, we shall study the following problem || 1C, 11 1: we are given a time series 
generated by F, and we want to estimate the maximal subgroup EgG under which 
the natural invariant measure (BRS-measure fll2"] ) is invariant, 

H(B) = n(hT x B) h G £ . (3) 

Here, B is any measurable set intersecting the attractor, and we assume that the 
initial conditions for the time series are //-random. Since this is a statistical decision 
problem, we cannot expect an unambiguous answer. In a proper treatment we should 
work out the detailed statistics of any used observable, and obtain the corresponding 
significance criteria. Following P, [10], |TTJ and a venerable tradition established in 
physics, we shall not do so here but hope that the results will be so clear cut that 
such formal significance criteria will not be necessary. 

In order to find the symmetry of a measure [i from a set of N points x^, k — 
1 . . . N drawn randomly from this measure, there are several possible strategies. 

(a) Visual inspection. For low dimensions, this is of course the method of choice. 
At higher dimensions, we can still visualize by making projections or cross sections. 
But usually, these will be less informative and more cumbersome. 

(b) We can partition phase space in some way such that each g G G acts as a 
permutation P g on the elements lAi of the partition, 

gU t =U 3 , j = P g i. (4) 

If the measure is symmetric, then we have 

lM(Ui) = »{U Pgi ). (5) 

Therefore, if the time series points are iid (identically and independently distributed), 
then 

rii « n Pgi ± 0(y/nl), (6) 

where n« is the number of points in U{. In simple cases this can be very efficient, 
but it becomes very unwieldy in high dimensions and if G-invariant partitions are 
tricky to construct. 

(c) We can consider averages over functions which themselves are not invariant 
under G. Assume that 0(x) is such a function. Then 

N 

(0) N' 1 £ 0(x fc ) (7) 
fc=i 

will in general be different from 

N 

(g<P) = N- 1 Y / <P(g- 1 Xk) (8) 
fe=i 
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unless [x is invariant under g. Functions 4> for which (0) = (g<p) only if fi is invariant 



are called detectives in ||, |l(| |TTJ (see these references for a more precise definition), 
and the numerical study of symmetry properties in these papers is based mainly on 
this method and on visual inspection of projections. 

(d) In the present paper, we shall use yet another method based on an idea 



by H. Kantz fl3| . In that work, the author was concerned with the more general 
problem of verifying numerically whether two measures are identical, making use 
of two random samples of points drawn from these measures. It uses the cross 
correlation sum, a generalization of the correlation sum used e.g. in estimates of the 
correlation dimension. 

Assume that the two sets of points are X = {x&, k = 1 . . . iV} and Y = {y^, k = 
1...N}. They are randomly drawn from the same space, but according to two 
different measures /i and v. The correlation sum for the set X is defined as 

where is the Heavyside step function. Thus Cxx measures the fraction of pairs 
having a distance < e. The cross correlation sum is defined similarly by 

^) = ^EEe( e -|x-y|). (io) 

Obviously, Cxy ~ Cxx up to statistical fluctuations if X and Y are distributed 
according to the same measure. Otherwise, we expect an inequality. More precisely, 
let us define the cross correlation integral by 

C^(e) = J d/*(x) J du(y) 6(6 - |x - y|) . (11) 

Almost surely, Cxy(e) — > C^ u (e) for iV — > oo, and Cxy(e) is an unbiased estimator 
for C M 2,(e) in the sense that C^ u (e) = (Cxy(e)) where the brackets denote the average 
over different random samples X and Y. For sufficiently small e and for absolutely 
continuous measures /i and v one can then show [T! 



ClM < C m (e) C w (e) . (12) 



Unfortunately, the same is not necessarily true for Cxy(e) with finite iV [13], due to 
the omission of "diagonal" terms in the double sum in eq.(|]). This omission is nec- 
essary since otherwise Cxx(e) would be a biased estimator for C w (e). Nevertheless, 
the analogon of eq.(^) for finite iV will be violated only in exceptional cases due to 
statistical fluctuations, and only if \i w v. 

Also, eq . (|12"D cannot be shown for all e and for arbitrary measures. This follows 



from Bochner's theorem and the fact that the Fourier transform of 0(e — |x|) is not 
positive definite. To avoid this problem, we could e.g. replace O by a Gaussian, 
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so that C^ v {t) = J dfi(x) J dv(y) e~^ x ' y ^ I 6 and similarly for the discrete cross 
correlation sums. We did not do this since this would make the numerics less efficient, 
and since we did not encounter any practical problems when using eq.flTOD in the 
applications described in sec. 3. 

Thus we use as a similarity measure the cross correlation ratio 

r(e) = - ° Me) . (13) 
^C xx (e)C YY (e) 

For very large e it is equal to 1 since Cxy(e) = Cxx{e) = Cyy(e) = 1 for sufficiently 
large e (we assume that the measures have compact support). On the other extreme, 
r(e) — > for e — > if the measures // and v have disjoint supports. The latter is 
typically || but not always || true if /i is a (time) invariant measure which is not 
invariant under some g G G, and v is the g-image of /i. For intermediate values 
of e we expect a monotonic increase, though we cannot exclude a more complex 
behavior. We can define a characteristic distance e c between the sets X and Y by 
demanding for instance r(e c ) = 1/2. 

In the present application, X is a time series and Y is just the image of X under 
a group action g, 

Y = gX = {<?x| x G X} . (14) 

Furthermore, we assume that the norm |.| is invariant under g (in the following, we 
use Euclidean norm). Thus Cyy = Cxx, arid the cross correlation ratio for element 
g is equal to 

r M = r h ■ ( 15 ) 

^xx{£) 

A problem arises from the fact that points in a time series are correlated and 
do thus not form a sample of iid points. This has mainly two more or less disjoint 
effects. The first is that it might take rather long until an orbit visits all parts of the 
attractor. This is particularly important in the vicinity of a bifurcation, and when 
there is intermittency or transients with very long characteristic times r. In such 
cases the only remedy is to take time sequences which span a total time T » r. 
To avoid huge data sets, one should then decimate the data. We thus take 

x n = x(nm dt) , n — 1, . . . N (16) 

where dt is the integration step and m is a sufficiently large integer. Independent 
from this, data on short scales might be correlated on time scales t < t corr with 
t corr > m dt. In this case, even the decimated data are correlated. As pointed out 



by Theiler |TJ]], it is then necessary to omit pairs in eqs. (p|JT0|) whose indices differ 
by less then n = t corr /m dt. Thus our final formula for r g (e) is 



TT N 


,9(6- 


|Xj — 9 X k\) @|i-fc|>n 




=l©(e- 


| x i X fc|) © \i— fc|>no 



r 9 (e) = Z . (17) 
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Finally let us discuss some immediate consequences following from the facts that 
Cxy is symmetric under the exchange X <-> Y, that G is a group, and that the 
norm is invariant under G. The first is that Cx, g x = CgX,x = Cx, g -^x for all g G G, 
and hence 

r g (e) = r g -i(e), 9 6G. (18) 

Next assume that h is in the symmetry group £ of the attractor. Then /iX « X 
and thus Cx, g x ~ Cx, g hx ~ Cx,h g x with errors vanishing for X — ► oo. This implies 
that r 9 takes the same value in the entire left and right cosets of £, 

= r ah(t) = r hg {e) , h e £ , g E G , (19) 
up to statistical fluctuations. 



2 Numerical Example 



Since we want to test our method against the method of detectives || [K], jlTf , we 
use exactly the same model as in [11]. It models the gait of a hexapodal animal, and 
is described in detail in 



It involves six identical nonlinear oscillators, coupled 
in a ring-like fashion. More precisely, the equations of motion are (note the misprint 
in the second line in fTTH ) 



dt 

dy P 
dt 



-Ax p + y p + (x 2 p + y 2 p ){Px p - Qy p )- 
A(4(x p _i + x p+x ) - 2(?/p_ 1 + Vp+x)) 



(20) 



-x 



p - 4% + (xp + yl){Py P + Q 



A(2(x p _x + x p+1 ) + 4(y p _! + y p+ i)) , p = 0, 1, ... 5 (mod 6) 



The system is equivariant under cyclic or anticyclic permutations of the indices of 
the oscillators. This would give the dihedral group D 6 as equivariance group. It 
consists of 12 permutations and is a subgroup of S12. Its generators are 



/v 



1 2 3 4 5 
5 4 3 2 1 



(21) 



C 



1 2 3 4 5 

1 2 3 4 5 



which was missed in |TT|. Thus the full equivariance group is 



G = D fi x Z 2 



(22) 



But in addition, the above system has a Z 2 symmetry (x p , y p ) — > (— x p , —y p ) 



(23) 



Attractors can have as symmetry groups any subgroup of G. Disregarding the 
factor Z 2 for the moment, the subgroups of D 6 are (up to conjugacies) 1, Z 2 = 
{1, C 3 }, Z 3 , Z 6 , D x = {1, k}, Eh = {1, <}, D 2 , D 3 , D 3 and D 6 . 
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In []IT| |, the attractor was projected for visual inspection onto a two dimensional 
subspace by means of the variables 



u = X\ — x 2 + x 4 — x 5 (24) 

and 

f = ~7=(^i + ^2 - 2x 3 + x 4 + x 5 - 2x ) . (25) 

In the (u, v) plane, the action of £ is a clockwise (negative) rotation by an angle 
27r/3, while k corresponds to a reflection on the axis u = 0. Notice that this is not a 
faithful representation of D 6 since ( 3 is represented by the unity (u, v) — > (u, v). The 
figures shown in |TTJ suggest that all attractors are symmetric under the additional 



Z2 symmetry. This was also verified by ourselves in all cases, whence we shall 
disregard this Z2 in the following and consider only symmetry under the elements 
ofD 6 . 

As distance measure in the correlation sums we used r = [J2 p (x p — x' p ) 2 + J2 p (y P — 

y' P ?] 1 ' 2 - 

Since our goal was to compare the correlation sum method to the method of 
detectives, we followed [[! ij in keeping the parameters P and A fixed at P = —5 and 
A = 1.2. The parameter Q was used as control parameter and was varied in the 



range 25.0 < Q < 27.2 studied also in [11 



We integrated eq.(]2"T|) with a fourth order Runge-Kutta algorithm with time 



step dt = 0.01. Again this was done in order to compare with |TT|. This time step 
might be too large for some of the Q- values considered, whence some of the observed 
structures might be artefacts of the time discretization. But this is not much of a 
concern: the time discretized system is as good a nonlinear dynamical system as the 
ODE, and studying its symmetries should be neither more easy nor more difficult. 

An important point where we departed from [|nj was the global time scales. We 
found that several of the attractors described in JTTJ] represented actually transients 



(e.g. for Q = 26.9, the final attractor is as shown in fig. 3 below, and not as described 
in |TT||). Thus we discarded in all transient of length t > 10000. In addition, 

we have the problem of long residence times on parts of the entire attractor discussed 
already in ||11|| . This gives rise to fake attractor asymmetries if the length of the time 



series is less than these residence times. To avoid this, we took only every 550th 
pointQ, i.e. m dt ~ 5 (see eq.([L6|)). Finally, we used a Theiler cutoff n = 100. 

Another point where we departed from [|TT|] was the number of points used. While 



the analysis of [] 1] was based on 10,000 - 50,000 points for each value of Q, we used 



only 5,000 points. This was done not only for minimizing the computational effort 



^or quasiperiodic attractors we found that large spurious correlations with very slow temporal 
decay could arise from near commensurability between the sampling frequency and some internal 
frequency of the attractor. To avoid this, we actually used an irregular sampling where we selected 
m randomly between 100 and 1000, such that (to) = 550. Without this trick, we would have 
obtained r g (e) >> 1 in several cases. 
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(which is higher for our method), but also to demonstrate that our method is very 
efficient already for small data samples. 

A last remark: as one should expect in general, we found in some cases that the 
system possessed several attractors which are not related by symmetry. For instance, 
depending on the initial conditions we found for Q = 25.7 either the attractor shown 
PHI or another one which closely resembles the one found in 111 at Q = 25.6. 



in 



3 Results 

In figures 1 to 4 we show four typical attractors together with the corresponding 
(cross-) correlation sums. For visualizing the attractors, we use the projection given 
by eqs.(|24]),(|25l). In each correlation sum plot, the number of pairs closer than e are 
plotted against e for each of the ten group elements g = 1, (, ( 2 , ( 3 , k, k(, k( 2 , k( 3 , 
k( a and k( 5 . Results for the remaining elements ( A and ( 5 are not shown since they 
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Figure 1: 



Panel (a): attractor for Q = 25.04 projected onto the (u, u)-plane by means 
Panel (b): log- log plot of the cross and auto correlation sums (without 
for Q = 25.04. The different symbols have the following 



of eqs.(|4],|25D 
normalization factors N 



meaning: O: g = 
g = K - * : g = K (. 



1 or k£ , +: g = ( or kQ 



3. 



□ : 



C or kC 



C or kQ°; A: 



coincide due to eq.(^) identically with those for ( 2 and (, respectively. Also, as we 
said already in the last section, we do not show the cross correlations for the inversion 
(x, y) — > (—x, — y) since all attractors were found to be symmetric under it. We show 
the numbers of pairs instead of the correlation sums with conventional normalization 
as in eq.(|9|), since the square roots of these numbers give rough estimates of statistical 
errors for small values of e where these errors are most important (for rigorous error 
estimates of autocorrelation sums, 



sec 



If an attractor is symmetric with respect to a group element g, we expect the 
corresponding Cx, g x to coincide with Cx,x- Within the expected statistical errors, 
the symmetries found in this way were also verified visually in the projections, except 
for g = ( 3 which is mapped by the projection (E3),(E3) onto the unit element. Thus 
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Figure 2: Same as fig.l, but for Q = 27.2. 



invariance under ( 3 cannot be read off from figures la to 4a, and can be seen only 
from the corresponding correlation sum plots. 





epsilon 



Figure 3: Same as fig.l, but for Q = 26.9. 



In fig.l, the projection suggests a symmetry under the inversion (u, v) — > (—u, —v), 
and maybe even under a rotation by 60°. We might be tempted to interpret this 
as invariance of the original attractor under Z 2 or Z 6 . From the correlation sums 
we see that this is not correct. Apart from the inversion (x, y) — > (— x, — y), the 
attractor has no symmetry at all, since none of the cross correlations agrees with 
the auto correlation sum (we notice that the invariance under (u, v) — ► (—it, —v) 
seen in all figures in JXTJ] is due to this inversion, and is not due to the particular 
projection used). 

The other extreme of a completely symmetric attractor is shown in fig. 2. We see 
that all cross correlation sums coincide. 

An example of D x invariance is shown in fig. 3. In the correlation sum plot we 
see four different curves. The one on top contains C x>x and Cx, K (x- The other 
curves contain C x ,(x, C x>kX , C x , K ^x (second from top), C Xt p x > C XM s x , C XM s x 
(third from top), and C X) ^ X , C XyKC _i X (bottom/rightmost curve). This degeneracy 
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epsilon 



Figure 4: Same as fig.l, but for Q = 26.1. 



of curves follows immediately from eq.flHj) and the detailed group structure of D 6 . 

Finally, we show in fig.4 a case with Z2 symmetry. We seem to see there three 
curves, but the rightmost one indeed consists of a triplet of lines which can only be 
resolved at higher resolution. Thus we have five curves altogether, each correspond- 
ing to two cross correlation sums. This is exactly as expected from eq. (|T9|) . The 
near degeneracy of the rightmost three lines is explained by the very weak breaking 
of D2 symmetry which is also seen in the fact that the two top curves are not very 
different except for very small e. 

In all examples (including many more which are not shown here), the symmetry 
suggested by the correlation sums is compatible with that seen by looking at the 
(u, v )-projection of the attractor. 

For a systematic investigation we have repeated this analysis for all values of 
Q e [25.0, 27.2] in steps of AQ = 0.01 (221 values altogether) . We repeated this 
three times with different types of initial conditions, in order to test for multiple 
attractors. In the first set of runs, we used the same initial values for all Q. In 
the second, different random initial conditions were used for each Q. Finally, in the 
third set of runs we used the last configuration from the simulation with Q as initial 
condition for simulating at Q + AQ. In some cases (such as Q = 2.57, e.g.) this 
gave clear indications of multiple attractors which are not simply rotated images of 
each other. 

Plots of cross correlation sum ratios against Q are shown in figures 5-10. These 
plots show results from the third type of initial conditions^, but the plots obtained 
with the other two types of inital conditions were very similar. Each figure corre- 
sponds to one group element. In each figure, three curves are shown. They represent 
r g {e) at those three values of e where the unnormalized autocorrelation sum was 
equal to 1000, 3000 and 10,000, respectively. In this way we guaranteed that the 
statistical errors were roughly independent of Q, although the size of the attractor 



2 For simplifying the figures, each attractor was first brought by a conjugacy into a standard 

k£" , n = . . . 5, was obtained at either 

g = k or g 



orientation such that the maximal r g for reflextions g 
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1.2 




25 25.5 26 26.5 27 

Q 

Figure 5: Cross correlation coefficients r g (e) for g = (, plotted against Q. Three data 
sets are superimposed, for distances e at which 1000, 3000 and 10,000 e-close pairs were 
found. The corresponding symbols are O (1000), + (3000), and □ (10,000). 





11 





Figure 9: Same as fig. 5, but for g = k(. 
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50000 100000 150000 200000 250000 300000 350000 400000 450000 500000 

time 

Figure 11: Cross correlation coefficient for g = k, Q = 25.22 and e = 0.36 plotted against 
time. For this, a very long time series was cut into pieces each containing 5000 points, and 
correlation coefficients were computed for each piece independently. On the horizontal axis 
is plotted physical time. In contrast to the previous plots, the average sampling frequency 
is « 1. 
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Figure 12: Final symmetry attributions versus Q. Notice that symmetries D2,D3 and 
D3 never appear. 



depends on it considerably. 

We do not show plots for g = k( 3 , k,( 4 and k( 5 since we found that r g = 1 for 
these elements if and only if r K( *2 = 1 (within statistical errors, of course). Thus for 
all values of Q the symmetry group Eq either contains all the elements k( 2 , k( 3 , k( 4 
and k( 5 , or none of them. This implies that none of the attractors had D2, D3, 
or D3 as its symmetry group, since all these groups (including all their conjugates) 
contain some elements of the set {k( 2 , k( 3 , k( 4 , k( 5 } but not all of them. 

For most values of Q this gave unique symmetry attributions, except for very 
few points mostly near Q w 25.2. In this range, the orbit switches occasionally 
between different geometric objects (which are of course not attractors in the strict 



sense), with vastly different time scales. This had already been noticed in [11]. As 



an example we show in fig. 11 the correlation sum ratio for g — k, e — 0.36 and 
Q = 25.22 as a function of time. For this figure, a very long trajectory (t = 5 x 10 5 
and 5 x 10 5 points as well) was partitioned into 100 pieces, and correlation sums 
were computed for each piece seperately. We see how r K switches beween and 1, 
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with a typical time scale of 10 . 

If this time scale is comparable to the time span covered by the set X, the latter 
looks as if obtained from a superposition of two or more measures. The same is 
obtained if part of the time series comes from a very long lived repellor. It is not 
easy to distinguish between these two possibilities using correlation sums only. But 
the presence of either of them is easily recognized in figs. 5-10, since r g (e) is roughly 
independent of e except for very small and very large values of the latter. 

Our final results are shown in fig. 12 where we have indicated the invariance 
group for every Q. Multiple entries result from coexisting attractors and/or from 
the above ambiguity. 



When comparing these results with those of [11], we observe several differences: 



(a) We see much more structure simply because we use a finer sampling in Q 
(AQ = 0.01 instead of 0.05 resp. 0.1). For that reason some of the very nar- 
row windows as those at Q = 25.92 and at Q = 25.96 — 25.97 with £ = 1 and Z 2 
were simply missed in |[L1| . 

(b) The window at Q — 26.89 — 26.92 with symmetries 1 and JDi was presumably 
missed in |L1| since it becomes visible only after a very long transient has died out. 
A long transient (and not intermittency as suggested in |11| ) is also the reason for 
the transition from D 2 to T)q at Q = 26.15 found in []TT ] . 



(c) The same trivial reason might also explain why the large region with Z 3 sym- 
metry at 25.74 < Q < 25.84 has been missed in fllTf . But there is also the more 
interesting possibility that it was missed since the attractor there is nearly (but not 
quite!) symmetric under the "rotation" g = (. As we can see from fig. 5, we have 
to go to rather small distances to see unambiguously that (e) ^ 1, and this might 
easily be overlooked when using detectives. 

(d) For most points in the interval 25.18 < Q < 25.32 we find the intermittent 
behavior illutrated in fig. 11 which suggests that the actual symmetry is the full D 6 . 
In contrast, the authors of (Tl| quote Z$ for Q = 25.2 and Q = 25.3, which is the 
symmetry observed during the (short) times when the curve in fig. 11 is below 1 (see 
also the discussion in [11] concerning Q = 25.55). 

(e) For 25.35 < Q < 25.52 we observe a large window with very clean Z 2 symmetry, 
while [[ll[] quotes Z 6 also for Q = 25 A. The reason for this discrepancy is not clear 
to us. There is a Z 6 -symmetric transient, but its life time t tra ns ~ 100 should be too 
small to have created any confusion. 



4 Conclusion 

We have shown that cross correlation sums are very efficient for studying symmetries 
of strange attractors. They are numerically more cumbersome than the "detectives" 
used in [[U], ||, [10], [!J , but they can be used for smaller data sets and they give much 



more detailed information. In particular, they allow easily to estimate the amount 
by which a symmetry is violated. 
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We tested our method on a model with D 6 symmetry studied by means of detec- 



tives in |TT[ . We succeeded to obtain unique symmetry attributions even for control 



parameter values where the authors of ['1 1] failed in spite of much larger data sets. 
Part of the reason for this is due to rather trivial modifications which are indepen- 
dent of the method used for symmetry detection. These include discarding very 
long transient and using decimated time series (eventually with irregular time gaps, 
in order not to be influenced by commensurability effects between the sampling fre- 
quency and some periodicity of the attractor). But even taking this into account, it 
seems that our method makes more complete use of the information contained in a 
given time series. 

Just as the authors of [[ll|], we did not bother with formal statstical error es- 
timates. Nevertheless, and in spite of the small samples used (5000 points), we 
feel that purely statistical errors were not a serious problem in the present appli- 
cation. Whenever we found an ambiguity, it was due to systematic effects (long 
transients and/or huge intermittency time scales) which cannot be estimated reli- 
ably by any method. At least, our method issues a warning in such a case in form 
of near-constant cross correlation ratios. 

In addition, of course, the estimation of correlation sums allows also to study 
scaling properties of the attractor. Usually, correlation sums are computed only 
for that purpose, i.e. in order to estimate attractor dimensions. We have not yet 
mentioned them here on purpose, to stress that correlation sums contain much more 
useful information than just dimensions. Nevertheless, we might take the chance 
here to say that the attractor dimensions indicated by figs. 1-4 are between 2 and 
5. 
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